Phase equilibrium in two orbital model under magnetic field 
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The phase equilibrium in manganites under magnetic field is studied using a two orbital model, 
based on the equivalent chemical potential principle for the competitive phases. We focus on the 
magnetic field induced melting process of CE phase in half-doped manganites. It is predicted that 
the homogenous CE phase begins to decompose into coexisting ferromagnetic phase and CE phase 
once the magnetic field exceeds the threshold field. In a more quantitative way, the volume fractions 
of the two competitive phases in the phase separation regime are evaluated. 
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Manganites, a typical class of strong correlated elec- 
tron systems, have been intensively studied in the last 
decade, due to their unusual behaviors of potential appli- 
cations such as colossal magnetoresistance (CMR). The 
double exchange (DE) mechanism can explain the mag- 
netic transition qualitatively, but more complex mecha- 
nism responsible for CMR is not yet fully understood. 
The idea of phase separation (PS) was recently proposed 
to understand the physics essence underlying the amaz- 
ing behaviors of manganites, while more and more theo- 
retical and experimental evidences confirm the existence 
of PS due to the intrinsic inhomogeneity Q, 0, • 

Former investigations on the phase diagram of man- 
ganites revealed the first-order character of phase transi- 
tions between various phases, e.g. charge-ordered (CO) 
insulator and ferromagnetic (FM) metal 0- Q 0- E3| • The 
insulator-metal transition in manganites can be reason- 
ably understood as the consequence of percolation of FM 
metal filaments embedded in the insulated matrix, and 
there are plenty of experimental evidences to support this 
PS framework |2,|7|,|8j- Current theories on manganites 
mainly stem from the competition between several inter- 
actions: DE, super exchange, Hund coupling, electron- 
phonon interaction and Coulomb interaction |2j ■ Besides, 
the effect of quench disorder on PS dynamics is high- 
lighted, especially on the large scale PS. The theoretical 
progress has enabled us to sketch the phase diagram in 
some special regimes from calculation and identify the PS 
regime in parameter space with various microscopic mod- 
els [UEJEI]- Nevertheless, it is still unclear theoretically 
how the PS develops, especially under an external pertur- 
bation, e.g. magnetic or electric field. In other word, it is 
of interest to not only identify the existence of PS regime, 
but also concern how the phase separation occurs and how 
it evolves upon external perturbation, because potential 
applications call for more sufficient theoretical interpre- 
tation. For instance, the CMR effect, which is one of the 
most attracting topics in the physics of manganites, may 
be described by the resistor network model phenomeno- 
logically, based on the percolation mechanism. In such 
case the volume fraction of metal phase is the key input 



variable, which, however, lacks credible theoretical inves- 
tigation yet. In earlier studies, this important variable 
was obtained from experiment or toy model Il2l . 

In this letter, we attempt to study the phase equilib- 
rium (PE) in a two orbital model. We emphasize par- 
ticularly the evolvement of PS upon the magnetic field 
perturbation. So we call it PE instead of PS in this work. 
Let's begin with a simplified model Hamiltonian which 
has been frequently used Q: 

Hami. = - ^ ^7'^^ -Jfl^s/S, 

+j AF J2 Si-Sj + H-J^gisi + Si)^) 

<i,j> i 

where a. is the vector connecting nearest-neighbor (NN) 
sites and dl lcr (di +ai ',j) is the generation (annihilation) 
operator for e g electron with spin a in the 7(7')-orbital 
on site i(i + a); t" r , is the NN hopping amplitude be- 
tween 7 and 7'-orbital (d x 2_ y 2 as a orbital, d^ z i_ r -2 as 
b orbital) along a-direction, with t\ h — to > (to is 
taken as energy unit), t x aa = t\ a = |i , t x hh = t\ b = \t Q , 

4-V 4-V 4-x j-x \/3 -f j-z j-z j-z n 

L ab ~ l ba ~ L ab ~~ L ba ~ 4 L °> L aa ~ L ab ~ L ba " u ' 

respectively; Si is the spin operators for t 2g core on site 
i, while Si for e g itinerant electron. The first term repre- 
sents the kinetic energy (DE process) which leads to FM 
spin arrangement. The second term is the Hund coupling 
of e g and t 2g electrons where Jh > is large enough to 
be regarded as infinite, so the spin of e g electron is always 
parallel with the same-site t 2g spin. The super exchange 
interaction Jaf > prefers to couple NN t 2g spins anti- 
ferromagnetically; The last term represents the magnetic 
field contribution with magnetic field H and Lande factor 
g. The electron-phonon coupling and Coulomb repulsion 
are not taken into account in this Hamiltonian and their 
effect on PE will be discussed below. 

The above simplified Hamiltonian can be solved ex- 
actly once a prior t 2g spin pattern is given. In real man- 
ganites, various t 2g spin patterns exist corresponding to 
abundant phases, e.g. FM, antiferromagnetic (AFM), 
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FIG. 1: (color online) (a)Density of state of FM/CE/CAFM 
phase. Here only regime E < is displayed since the symmet- 
rical E > part is empty in ground state. In GAFM, there 
is only a single energy level at because of the localization 
of e 9 electron, (b) Chemical potential /j, as a function of e g 
electron concentration. The energy gap to for the CE phase 
at n — 0.5 is evident both in (a) and (b). 



CO, orbital-ordered (00) phases. In this work, sev- 
eral typical tig patterns confirmed from experiments: C- 
type AFM (CAFM), G-type AFM (GAFM), FM and 
CE phase are chosen as the candidates to compare. 
The CAFM phase is constructed by antiferromagnet- 
ically coupled one-dimension FM lines, while the CE 
phase is constructed by antiferromagnetically coupled 
one-dimension zigzag FM chains and is found to be 
CO/OO 0E[Ii|. The GAFM takes the familiar AFM 
arrangement in all three directions. Then the Hamilto- 
nian can be exactly solved when the Hund factor Jh 
is simplified as infinite. The procedure of derivation is 
straightforward and the details can be found in Ref.|2j. 
Then density of state (DOS, D{E)) of the these phases 
can be calculated from the dispersion relationship using 
numerical method, as shown in Fig^a). In addition, the 
chemical potential fj, of these phases is obtained simulta- 
neously by integrating the DOS, as shown in Fig^b). 
Consequently, the ground state energy is calculated. For 
instance, the energy of FM phase can be written as: 



E(n,H) _ 1 
N ~~ N 



27 

D(E)EdE + —Jaf - 3H (2) 



here N is the number of whole sites and is infinite ideally. 
The first integral term gives the energy of all e g electrons; 
term 27 Jaf I '2 arises from the six NN FM correlation 
between \S\ = 3/2 ti g cores in classical approximation; 
factor 3 before H is calculated by multiplying ti g spin 
3/2 with Lande factor 2. Besides, the influence of H 
on D(E) has also been taken into account. The average 
e g electron concentration n € [0, 1]. The DOS, chemical 
potential and ground energy of other phases can also be 
calculated exactly just as the FM case. 

However, what should be noted is that none of these 
phases can be stable over the whole doping range. In 
order to determine which phase is the preferred one at 
a given concentration, the ground state energy of these 
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FIG. 2: (color online) System energy E (ground state under 
zero field) as a function of e g electron concentration n, with 
Jaf ~ 0.0278to- Upon increasing e g electrons n, a transition 
sequence of the ground state: GAFM^CAFM^CE^FM, is 
shown. The CE phase as ground state is possible only in the 
narrow range around the half-filling point. 



phases, when Jaf is set as 0. 0278io and H = 0, is plot- 
ted in Fig|2 The preferred ground state is the energy 
minimal state. It should be GAFM as n ~ because 
the interaction is almost pure AFM super exchange. As 
n ~ 0.3, the preferred state is CAFM. The CE phase 
can appear only in a narrow regime n ~ 0.5. When the 
gain from kinetic energy suppresses the loss of super ex- 
change energy in the large n regime, the FM becomes the 
stable phase. With this calculation, the phase diagram 
over the whole concentration regime has be developed. 
Of course, the phase diagram of real manganites is more 
complex than this simple sketch for two reasons: first, 
the Hamiltonian Eq.(l) is oversimplified and secondly the 
candidate phases chosen here are not complete but only 
four phases. Even though, the calculated phase diagram 
is quite similar to that of some typical manganites (zero- 
temperature): e.g. Ndi-^Sr^MnOa (here n = l — x) |17j . 
In fact, it is shown that the phase diagram of other man- 
ganites: e.g. Lai-ajSr^MnOs or Pri_ x Ca x Mn03, can be 
reproduced roughly by adjusting the value of Jaf, whose 
role will be revisited below. An important truth revealed 
here is that no matter what Jaf, the CE phase can ei- 
ther appear in the narrow regime n ~ 0.5 or simply be 
unstable over the whole concentration regime. 

The above theoretical approach based on Hamiltonian 
Eq.(l), whose origin can be found from Ref. 2j, allows us 
to study the PE in PS regime. Here, a simple but repre- 
sentative case: the PE between FM phase and CE phase 
with n = 0.5, will be studied. It corresponds to the melt- 
ing process of a type of CO state (here it is CE phase) un- 
der external magnetic field. Since for FM metallic phase 
the lattice distortions are absent and the e g electrons are 
delocalized 0, the electron-phonon coupling and on- 
site Coulomb repulsion are unimportant. On the other 
hand, the CE phase can also be reasonably described as 
a band insulator by this Hamiltonian at n = 0.5 case 
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HE!]. Therefore, Eq.(l) is suitable to deal with the PE 
between FM metallic phase and CE phase, noting that 
Eq.(l) can be exactly solved without scale issue. 

The PE principle for a PS system can be represented by 
the equivalence in chemical potential between the com- 
petitive phases, i.e. FM and CE to be considered. From 
the above calculation (Fig^b)), it is seen that at n = 0.5 
and under zero field, the chemical potential of pure FM 
phase is about — 1.52£ and that of CE phase is — 1 . 
Therefore, the chemical potential [i for the possible PS 
system would be in the range [— 1.52i , ~t ]. The sum of 
e g electrons in the FM and CE phases can be calculated: 

nA = ~w I D{E)dE (3) 

here pa is the volume fraction of A (FM/CE) phase. The 
upper limit of the integral should meet the equivalent 
chemical potential condition: \Ifm — fJ.cE — Then 
the following set of equations yields: 

f PFM + PCE = 1 / 4 n 

\ n FM + n C E = 0.5 

By numerical method, Eq.(0} can be solved with the 
prior assumed fi. Given parameters Jaf — 0.0278io and 
H = 0, the calculated relative volume fraction of FM/CE 
phase as a function /i are shown in Fig^a) (right axis). 
Then by including the factor pa into the energy equation, 
the total system energy: E = pfmEfm + PceEce can 
been calculated, as shown in Fig|^a) (left axis). It is seen 
that both pfm and E are monotonously decreasing func- 
tions of \i. The ground state must be homogenous CE 
phase because the PS state goes against energy. In case 
of nonzero magnetic field (e.g. H = 0.0I5io), the above 
calculation is repeated by taking the magnetic field con- 
tribution to DOS and energy into account. The results 
are plotted in Fig|3fb). Quite interestingly, a nonzero 
field results in a minimal in the E — fj, curve (here at 
/i ~ — 1.13£o)- The PS state associated with this minimal 
energy point is more stable with respect to the homoge- 
nous FM or CE phase. The value — 1.13to is the real 
chemical potential of this PS ground state which consists 
of 57% CE phase and 43% FM phase. 

By varying the value of H, one can repeat the above 
calculation and then obtain the relative volume frac- 
tions of the two phases as a function of H, as shown 
in Fig^Ja). When H is low (below the lower thresh- 
old H min ~ 0.0069to), the CE phase is robust against 
magnetic field perturbation, indicating the stable and 
homogenous CE phase as the ground state. Upon an in- 
creasing of H beyond H min , part of the CE phase begins 
to melt into FM phase, and pce is deceasing with increas- 
ing H under the equivalent chemical potential condition. 
In details, Pfm increases rapidly once H > H rn i n , and 
the growth becomes slower when H is further higher, 
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FIG. 3: (color online) Volume fraction of FM/CE phase (right 
axis) and corresponding total energy system E (left axis) as 
a function of chemical potential fi. Jaf = 0.0278io- (a)Zero 
field case, where the ground state is CE phase without phase 
separation because the energy decreases monotonously with 
chemical potential; (b)under magnetic field H = 0.015to, 
where a energy minimal appears at fi = — 1.13to- The co- 
existence with the fraction of FM and CE phase in this point 
is the ground state. 

as shown in FigQJa). Note here that the threshold of 
percolation for FM phase can be easily surpassed under 
a field slightly higher than H m i n . For instance, only a 
field of H = 0.009io is needed to obtain a FM phase 
of 24.7% (threshold for a three-dimensional simple cu- 
bic bond percolation), beyond which an insulator- metal 
transition may be expected. When H is extremely high 
(not shown in FigQJa)), the equivalent chemical poten- 
tial condition can no longer be satisfied, indicating a ter- 
mination of the PS state and the ground state will be 
homogenous FM state. 

The above calculation indicates that the lower thresh- 
old of magnetic field to melt the CE phase is 0.0069£o- 
Considering that to in low-bandwidth manganites is 
small, e.g. O.ley, the calculated H min is about 12T, 
a value consistent with the experimental data 0, |2(], 
EH . Note here that H min is not equal the direct en- 
ergy gap between FM and CE phase, which is one order 
of magnitude larger than experimental value (shown in 
FigEl about 0.08<o)- It means that the required magnetic 
field to destroy the CO insulator is strongly reduced by 
PS which can occur once the energy of competitive phases 
is close to each other. On the other hand, as identified 
earlier, parameter Jaf plays a key role in PS although it 
is the least intrinsic interaction in manganites • Other 
than the above case of field-induced sequences, abundant 
phenomena associated with PS in manganites can be pre- 
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FIG. 4: (a)Melting process of the CE phase under magnetic 
field. When the field is below the threshold H m i„, the ho- 
mogenous CE phase is robust. Then phase separation occurs 
when H > H m i n . (b) Phase transitions induced by Jaf un- 
der zero-field condition: FM phase at weak Jaf', coexistence 
of FM and CE phase at middle Jaf', CE phase at strong Jaf- 



dieted by our model through adjusting Jaf- For exam- 
ple, again at n = 0.5, a coexistence of FM phase and CE 
phase under zero field is predicted at 0.019to < Jaf < 
0.026^0- When Jaf is further reduced, a homogeneous 
FM phase as the ground state is possible even under zero 
field. The transition between three regimes: homoge- 
nous FM phase to PS state to homogenous CE phase, are 
identified in Fig0Jb). These transitions can be mapped 
to real manganites of wide-band to those of middle-band 
and then narrow-band [2|. 

It should be mentioned that the only parameter ad- 
justable here is the super exchange Jaf, and the energy 
difference between different phases is dependent on ra- 
tio Jaf /to- This is obviously oversimplified, referring 
to real manganites materials in which not only the dou- 
ble/super exchange but also the Jahn- Teller distortion 
and Coulomb repulsion play important roles. For in- 
stance, the Jahn- Teller distortion in CE phase will af- 
fect the energy band and DOS I n addition, al- 
though the phase diagram given in Figgis quite similar 
to those for some manganites, there are still some blem- 
ishes. For instance, a prediction of the correct concen- 
tration n corresponding to the A-type AFM observed in 
LaMn0 3 (n ~ 1) Q or Ndi^Sr^MnOa (x ~ 0.55) [13 
can not be given by the present model. However, if a com- 
plete Hamiltonian is employed, the calculation has to be 
oversimplified, e.g. limited in a small cluster which is dis- 
advantageous to deal with PS. Fortunately, Eq.(l) in the 
present work can describe FM/CE phase to some satisfac- 
tory extent and it is a good starting point to investigate 
the PE issue in PS systems against external magnetic 
field perturbation. Furthermore, the present approach 
represents a general roadmap to investigate the PE is- 
sues in manganites: e.g. phase competition other than 
FM-CE, of more than two phases, of different e g concen- 
trations, and effect of other perturbation than magnetic 
field etc. The key condition is the equivalence of chemical 



potential between competitive phases. 

In summary, the principle of chemical potential equiv- 
alence has been introduced to investigate the phase equi- 
librium of half-doped manganites under external mag- 
netic field. By employing the two orbital model, we have 
presented an explicit solution to the phase equilibrium 
between FM phase and CE phase. The magnetic field 
threshold required for melting of the CE phase has been 
calculated, consistent with the experimental results. The 
volume fractions of the two competitive phases in the 
phase separation regime as a function of external mag- 
netic field have been evaluated. In addition, the super 
exchange modulated transitions between ferromagnetic, 
phase separated and CE states under zero-field, is pre- 
dicted. 
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